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Abstract 


Woven  mesh  structures,  consisting  of  bonded  laminates  of  two-dimensional  plain-weave 
conductive  screens,  or  three-dimensional  orthogonal  weaves  are  described.  A  mathematical 
model  of  the  thermal  performance  of  such  a  mesh,  deployed  as  a  heat  exchange  surface,  is 
developed.  The  model  shows  that  the  effective  conductance  of  a  thin,  porous  heat  exchange 
matrix  is  proportional  to  the  square  root  of  the  product  of  the  specific  surface  area,  effective 
thermal  conductivity  and  mesh  heat  transfer  coefficient.  Furthermore,  we  propose  that  mesh  heat 
exchange  matrices  be  laid  out  in  a  serpentine  pattern,  similar  to  the  construction  of  a  conventional 
automotive  air  filter.  This  configuration  maximizes  heat  transfer  surface  area  while  it  minimizes 
coolant  pressure  drop  (due  to  the  reduced  superficial  mass  velocity).  Analysis  of  a  screen-fin  heat 
sink  that  employs  this  technology  indicates  that  the  optimal  convergence/divergence  angle  of  the 
serpentine  layout  is  »  20°  -  30°. 

Geometric  equations  show  that  these  porous  matrices  can  be  fabricated  to  have  a  wide  range  of 
porosity  and  specific  surface  area,  with  the  specific  surface  area-filament  diameter  product  of  3D 
stacked  weaves  approaching  2n.  Furthermore,  a  highly  anisotropic  thermal  conductivity  vector 
can  be  achieved,  with  the  in-plane  effective  thermal  conductivity  ranging  to  78.5%  of  the  filament 
material  values. 

Measurements  of  pressure  drop  and  overall  heat  transfer  rate  are  reported  and  used  with  the 
performance  model  to  develop  correlation  equations  of  mesh  friction  factor  and  Colburn  j-factor 
as  a  function  of  coolant  properties,  mesh  characteristics  and  flow  rate  through  the  mesh. 
Experiments  with  air  and  water  show  that  mesh  heat  transfer  coefficients  and  friction  factors  are 
comparable  to  those  achieved  with  other  expanded  materials.  However,  high  p-values,  coupled 
with  high  effective  thermal  conductivity  result  in  exchange  matrices  that  out-perform  other 
exchange  matrix  configurations.  This  is  particularly  true  for  exchange  matrix  configurations  that 
exploit  the  higher  effective  thermal  conductivity  achievable  with  this  technology. 
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Introduction 


Kays  and  London  [1984]  have  pointed  out  that  a  most  effective  way  to  increase  the  performance 
of  a  heat  exchanger  is  to  increase  its  surface  area  to  volume  ratio,  p.  Small-particle  packed  beds 
and  foamed  metals  are  porous  materials  having  large  /^values.  Unfortunately,  due  to  the 
tortuosity  effect  in  conjunction  with  the  high  porosity  (s)  of  these  materials,  their  effective 
thermal  conductivity  (ke)  is  relatively  small  so  that  much  of  the  gain  in  performance  obtained  by 
having  a  large  P  is  lost  by  having  a  relatively  small  ke.  Typical  values  of  effective  thermal 
conductivity  in  fused-particle  packed  beds  are  10%  -  15%  of  the  particle  thermal  conductivity. 
Commercially  available  metal  foam  such  as  aluminum  foam  has  an  effective  thermal  conductivity 
that  ranges  from  only  2%  to  6%  of  the  base  metal  value  [Calmidi  and  Mahajan,  1999]. 

An  anisotropic  porous  matrix  having  a  large  specific  surface  area  and  a  large  effective  thermal 
conductivity  in  a  particular  direction  will  result  in  a  very  effective  heat  exchange  surface.  A 
woven  mesh  of  heat  conducting  filaments  can  be  configured  to  have  these  characteristics. 
Geometric  equations  show  that  these  porous  matrices  can  be  fabricated  to  have  a  wide  range  of 
porosity  and  specific  surface  area;  and  a  highly  anisotropic  thermal  conductivity  vector  can  be 
achieved. 

These  attributes  allow  for  the  design  of  small,  high-performance  single-fluid  parallel-plate  heat 
exchangers  that  are  more  universally  applicable  than  conventional  heat  exchangers  because  the 
mesh  can  readily  be  made  to  conform  to  complex  surface  contours.  Because  of  the  high  thermal 
conductance  achieved,  exchangers  can  be  designed  for  applications  where  spatial  temperature 
uniformity  and/or  high  localized  spot  cooling  are  required. 

The  present  work  focuses  on  two  woven  mesh  configurations: 

1 .  Laminates  of  plain-weave  metallic  screens  stacked  to  form  a  three  dimensional  structure, 
and 

2.  Three-dimensional  orthogonal  weaves  that  consist  of  three  metallic  wire  filaments  having 
their  axes  aligned  with  the  coordinate  axes. 

The  fabrication  methodology  for  each  structure  is  explored.  Models  of  the  porosity,  specific 
surface  area,  effective  thermal  conductivity  and  overall  thermal  performance  of  the  mesh, 
deployed  as  a  single-fluid  heat  exchanger,  are  developed.  Measurements  are  reported  for 
fabricated  test  articles;  and,  mesh  Stanton  number  and  friction  factor  correlations  are  reported. 
The  thermal  performance  of  the  mesh  configurations,  deployed  as  heat  exchanger  surfaces,  are 
evaluated  in  the  light  of  other  surfaces.  Prototype  heat  exchangers  are  constructed  and  tested. 
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Fig.  1  Heat  exchange  implementation  of  a 
thermally  conductive  woven  mesh. 


Heat  Exchanger  Implementation 

Figure  1  shows  a  woven  mesh  implemented  as 

a  heat  transfer  surface  in  a  parallel  plate 
exchanger.  The  mesh  is  shown  in  edge-view 
in  a  channel  having  half-height  H.  The  mesh 
thickness  is  t.  Heat  (q)  is  conducted  from  the 
heated  plates  (at  temperature  Tb)  into  the 
mesh,  and  then  by  convection  to  the  fluid 
flowing  through  the  mesh.  The  mesh  wire 
filament  axes  are  arranged  with  the 
predominant  y-filaments  perpendicular  to  the 
channel  walls  so  that  conduction  from  the 
walls  is  facilitated.  Other  wire  filaments  in 
thermal  contact  with  the  predominant 
filaments  act  as  fins,  increasing  the  specific 
surface  area  of  the  mesh.  G{  and  7#  are  the  coolant  approach  mass  velocity  and  temperature, 
respectively.  AP  is  the  pressure  drop  across  the  exchange  matrix.  The  right-hand  segment  of  the 
figure  shows  the  expected  temperature  distribution  of  the  woven  mesh  (Ts)  and  the  exit-plane 
coolant  (Tf0).  2<f> 

Figure  2  shows  a  plan  view  of  a  segment  of 
the  coolant  passage  of  a  plate-fin  heat 
exchanger.  Fin  stock  normally  found  between 
the  plates  is  replaced  with  a  woven-mesh 
porous  wall  that  is  laid  out  in  a  serpentine 
pattern.  Mesh  wire  filaments  are  oriented 
between  the  side  plates  of  the  exchanger  as  in 
Fig.  1.  The  woven-mesh  porous  wall  segment 
so  formed  is  shown  in  plan  view,  and  typical 
flow  paths  for  some  fluid  elements  are  shown. 

This  configuration  results  in  a  relatively  low 
coolant  pressure  drop  since  fluid  parcel  paths 
through  the  mesh  are  short,  and  the  serpentine 

mesh  layout  results  in  a  large  face  area  leading  to  low  superficial  mass  velocities  of  coolant 
through  the  mesh.  The  divergence/convergence  angle  of  wall  segments,  2  <f>  controls  the  flow 


Porous 

Mesh 


Coolant  Flow 

Fig.  2  Plan  view  of  woven  mesh  serpentine 
wall.  Divergence  angle  controls  uniformity 
of  flow  through  wall  segments. 
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uniformity  through  each  segment.  Heat  transfer  is  by  conduction  through  the  module  side  plates, 
and  then  by  conduction  and  convection  within  the  (porous)  woven  mesh. 


Thermal  Performance  Model 

The  mesh  acts  as  a  porous  wall  of  porosity,  e  and  thickness,  t.  The  heat  transfer  rate  is  given  by 

q  =  U(tW)\rh-Tf,\  (D 


where  U  is  the  effective  conductance  of  the  porous  wall  and  (tW)  is  the  porous  wall  base  area. 

The  porous  wall  effective  conductance,  U  can  be  related  to  the  thermal  and  physical 
characteristics  of  the  woven  structure.  The  fluid  flow  path  length  through  the  porous  wall  is  short 
and  flow  rates  are  relatively  high,  so  local  thermal  equilibrium  between  the  fluid  and  solid  phases 
is  probably  not  achieved;  a  two  energy  equation  model  is  called  for.  Wirtz  [1997]  assumes  that 
the  solid  phase  temperature  is  only  a  function  of  y,  and  the  local  heat  flux  between  the  fluid  and 
solid  phases  is  characterized  by  Newton’s  cooling  law 

=  (2) 


where  h  is  the  mesh  heat  transfer  coefficient.  Eq.  (2)  couples  the  solid  and  fluid  phase  energy 
equations. 

An  energy  balance  on  a  dy- 1  slice  of  solid  phase  material,  which  balances  the  increase  in  internal 
energy  of  the  fluid  phase  with  the  net  conduction  in  the  solid  phase,  leads  to  the  following  solid- 
phase  energy  equation 


d% 

dy 2 


cGi 

keyt 


VfAy)-TfA 


(3) 


dT. 


with  boundary  conditions:  Ts(o)  =  Tb  and  — 

dy 


=  0  (symmetrical  heating  assumed).  An  energy 


\y=H 


balance  on  a  dx-  dy  element  of  fluid,  which  introduces  the  particle  heat  transfer  coefficient 
through  eq.  (2),  gives 


dr  pSt  n 
dx  £ 


(d) 


sh 


where  r  =  T  (y)  - Tf  (x, y) ,  with  boundary  condition  r(0,y)  =  T  (y) -Tn.  St  =  — —  is  the 

cG' 

mesh  Stanton  number.  Wirtz  [1997]  assumes  that  Tf(x,y)  =  (7),  +  Tf  o)l2 .  This  assumption  leads 
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to  closure  of  the  two  energy  equations.  The  following  is  an  extension  of  the  above-described 
model,  which  does  not  make  this  assumption. 

Assume  the  porous  wall  has  uniform  porosity  and  thickness.  Then,  the  fluid  flowing  through  the 
porous  wall  will  follow  parallel  streamlines  (y  =  constant)  so  that  a  parcel  of  fluid  will  be  exposed 
to  solid  phase  material  at  a  fixed  temperature.  Further  assume  that  h  =  constant.  Then  eq.  (4)  can 
be  integrated  across  the  porous  wall  to  give 

Tf.o(y)-TfJ  =[7;O0— (5) 

with  ntu  =  St the  number  of  transfer  units  of  the  mesh.  Substitution  of  eq.  (5)  into  eq.  (3) 
e 

and  integration  gives  the  solution  for  the  solid  phase  temperature  distribution 


Ts (y)  -Tf  i  _  cosh[w(// - y)] 

Tb  -  Tf  i  cosh[w//] 

(6) 

i 

i 

1 

^1— 

■^r 

H 

S 

(7) 

where  G,.  = 


cG,#2 

keyt 


is  a  dimensionless  coolant  superficial  mass  velocity.  Then,  the  heat  flux 


dT 

ql  =  -he — -  may  be  determined  by  differentiating  eq.  (6).  Substitution  of  the  result  into  eq. 

toy.  o' 

(7)  gives 


—  =  /«//•  tanh(m//) 
keu 


(8) 


Finally,  recognizing  that  the  maximum  heat  transfer  rate  occurs  when  all  the  coolant  is  heated  to 
the  base  temperature  of  the  porous  wall,  we  can  write  an  expression  for  the  effectiveness  of  the 
screen  laminate  fin 


4*  (max) 


tanh(/w//) 


(9) 


Figure  3  compares  the  effectiveness  of  the  present  model,  eq.  (9)  with  that  of  Wirtz  [1997]. 
The  figure  plots  tihe  exchange  matrix  effectiveness  versus  number  of  transfer  units  for  two 
dimensionless  superficial  mass  velocities.  The  figure  shows  that  the  new  model  gives  a  smooth 
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variation  in  effectiveness  as  a  function  of  ntu, 
and  it  predicts  rj  somewhat  lower  than  the 
previous  model.  The  effectiveness  is  seen  to 
asymptote  to  unity  as  ntu  increases  when 

Gt->  0. 

Equations  (7)  and  (8)  show  that 

U  «  Jfikeyh  •  F{ntu)  (10) 

where,  for  the  current  application, 
F(ntu )  «  0(1) .  Therefore,  thermal 


ntu 

Fie.  3  Comnarison  of  norons  models. 


performance  is  maximized  by  increasing  the  specific  surface  area  (p),  the  effective  thermal 


conductivity  (key),  and  the  mesh  heat  transfer  coefficient  (  h  ). 
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Plain-Weave  Screen  Laminates 


Screen-Laminate  Geometry 
Figure  4  shows  two-dimensional  plain-weave 

screens  stacked  together  to  form  a  screen 

laminate.  Each  screen  has  woven  wires  of 

diameter  dy  and  dz,  with  axis  parallel  to  the  y- 

and  z-axis,  respectively.  Wire  spacing  is 

designated  by  the  mesh  numbers,  My  and  Mz. 

The  laminate  has  thickness, 

t  =  cf  •  n{dy  +  dz),  where  n  is  the  number  of 

screen  layers  of  the  lamination,  and  cf  is  the  My 

compression  factor.  cf  accounts  for 

interleaving  of  wire  filaments  of  adjacent 

screens,  and  wire  crimping  at  wire 

.  ,  , .  ....  .  ,  j  .  ,  Fig.  4  Screen  laminate  geometry, 

intersections.  Wire  filaments  are  bonded  at  °  °  J 

intersections  to  facilitate  conduction;  and,  successive  screen  layers  are  bonded  together  to  add 

rigidity  to  the  structure.  Successive  screens  can  be  arranged  “ in-line ”  as  shown  in  Figure  4,  or  in 

staggered 5  configuration,  where  alternate  screen  layers  are  offset  in  they  and  z  directions  by  0.5 

M  ~J.  Screen  laminates  having  dy  =  dz  and  My  =  Mz  are  designated  isotropic.  If  the  angle 

between  wire  filament  axes  (the  weave  angle)  is  90°,  as  shown  in  Fig.  4,  the  screen  is  an 

orthogonal  weave.  Non-orthogonal  weaves  (diamond  weaves)  are  also  possible. 


Fabrication  Methodology 

Screen  laminate  fabrication  is  very  straightforward.  Commercial  grade  screen  is  stacked  to  form 
a  lamination  with  the  screen  elements  either  aligned  or  staggered.  The  wire  filaments  are  then 


dip-brazed  or  re-flow  soldered  to  bond  the 
wire  filaments  at  their  intersections,  and  bond 
adjacent  screens  together.  The  resulting 
structure  is  veiy  stiff.  In  fact,  there  is  a 
commercial  stainless  steel  filter  product  called 
Rigimesh®. 

Figure  5  shows  multi-ply  prototypes  with  d  - 
0.432  mm  (0.017”),  M  =  20  in'1.  The  two 
different  screen  alignments  are  shown:  the  in- 


Fig.  5  Inline  and  staggered  plain-weave 
screen  laminates. 
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line  configuration  has  successive  screens  with  wire  filaments  aligned,  and  the  staggered 
configuration  has  screens  stacked  somewhat  randomly,  while  wire  filaments  are  aligned  with  the 
y-and  z-axis.  Commercial  grade  copper  screen  (ks  -  400  W/mK)  is  pre-coated  with  a  solder 
bearing  flux  (95%  Sn/5%  Pb,  ks  =  5 5-60 W/mK).  Screen  layers  are  stacked  in  a  fixture,  and  the 
assembly  is  re-flow  soldered  at  230°C.  The  inline  configuration  has  s=  0.7,  /?=  2700  m"1,  and  hey 
-  50  W/mK.  The  staggered  configuration  has  s—  0.63,  j3  =  3350  m'1,  and  key  =  62  W/mK. 


Orthogonal-Weave  Porosity,  Specific  Surface  Area  and  Effective  Thermal  Conductivity 

Figure  6  shows  plan  and  edge  views  of  a 


section  of  screen.  Serpentine  wire  filaments 
have  diameter  dx  and  dy,  and  corresponding 
mesh  numbers  Mx  and  My.  The  wire  filament 

pitch  in  the  x-  and  y-directions  are  Mx~'  and 
My~l ,  respectively.  In  the  absence  of 
crimping,  a  single  screen  has  thickness  t/  =  dx 
+  dy.  Filament  lengths  are,  [Luo  and  Mitra, 


1999] 


Fig.6  Orthogonal  weave  unit  cell 


Fig.  7  Screen  stacking  configurations. 


Figure  7  shows  two  possible  arrangements  of  stacked  plain-weave  screens.  The  leftr  view  shows 
the  situation  where  successive  wire  filaments  are  aligned  so  that  the  thickness  of  a  laminate 
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consisting  of  n  screens  is  t„  =  n(dx + dy) .  The 
right  view  shows  the  situation  where 
successive  wire  filaments  are  not  necessarily 
in  line  so  there  is  some  interleaving  of  wire 
filaments  when  the  screens  are  stacked.  In  this 
case,  the  thickness  of  the  laminate  is 
/„  -n-Cf  (dx+dy)  where  c/is  the  compression 

factor.  The  magnitude  of  the  compression 
factor  depends  on  the  weave  pattern,  mesh 
numbers  and  wire  filament  diameters  of  the 
screen,  and  the  stacking  arrangement  of  the 
screen  layers.  If  we  restrict  our  attention  to 
screens  with  tx=  dx+dy ,  then  cf  <  1 . 


Fig.  8  Orthogonal  weave  transformed  unit 
cell. 


If  we  introduce  the  mesh  number  and  wire  filament  diameter  ratios  M  =  My!  Mx  and  d  =  dx/dy , 
then  it  can  be  shown  [Xu,  2001]  that  the  reduced  metal  fraction  is  a  function  of  [Mxdy,M,d)  and 
0  is  a  function  of  ((1  -  e),cf,M,d ). 

Koh  and  Fortini  [1973]  and  Chang  [1990]  have  reported  models  and  empirical  correlations 
specifically  directed  at  the  cross-plane  effective  thermal  conductivity  component  of  screen 
material.  We  have  found  no  reference  in  the  technical  literature  for  the  in-plane  component  of  the 
effective  thermal  conductivity. 
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5 

to 


Following  Chang  [1990],  we  transform  the 
segment  of  screen  unit  cell  shown  in  Fig.  6  to 
the  layered,  rectangular  cross  section  segment 
shown  in  the  figure  8.  Each  rectangular  wire 
filament  has  thickness  Tid/Ag  and  width  gd  so 
that  the  cross  section  area  of  each  filament  is 
7td2 1 A .  The  geometric  index,  g  is  a  measure 
of  the  contact  between  wire  filaments  at 
intersections. 

The  in-plane  effective  thermal  conductivity  in 
the  y-direction  may  be  determined  by 
considering  the  thermal  circuit  for  conduction 
across  the  transformed  unit  cell,  shown  in  Fig. 

9.  The  figure  shows  three  heat-flow  paths:  RS2 
is  the  thermal  resistance  for  conduction  along  the  axis  of  the  y-wire  filaments;  Rsi  is  the  thermal 
resistance  across  x-wire  filaments;  R/i  is  the  thermal  resistance  to  conduction  in  the  y-direction 
across  the  intervening  “fluid”  lying  between  x-wire  filaments;  and,  Rp  is  the  thermal  resistance  of 
the  fluid  slab  that  lays  between  the  y-wire  filaments.  These  resistances  are  given  as 


Fig.  9  Orthogonal  weave  thermal  circuit. 


R  -dx 

^  x-Ks-dx-Sx 


(16) 


Rs2  ~ 


8  Sy 


n-Ks-dy 2 


(17) 


Rfl  = 


4  g-Sy{My,dx)-4g2dx 
jc-Kt  • dx-Sx 


(18) 


4  g-Sy 

7T -Kf  \Sx-g-dy)-dy 


(19) 


We  further  require  that  the  volume  of  wire  filaments  be  preserved  across  the  geometric 
transformation  shown  in  Fig.  8.  Then  g  =  n  /(4  •  cf) .  The  Y-direction  conductance  across  the 
transformed  unit  cell  is. 


cf  ■ Sx-(dx+dy ) 
Sy 


(20) 
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where  R,ot  is  the  overall  thermal  resistance  of  the  thermal  circuit  (Fig.  9)  and  key  is  the  Y- 
direction  in-plane  effective  thermal  conductivity  of  the  screen.  Rearranging  Eq.  (20)  gives 


7  _mxMxdy(Kfs-\)  +  Ccf-Kfs  C-cf2-d-Kfs  /olx 

CJ’Key - r= — - +  7= r-r - —  — - = - v 

\d  +  \)D  \d  +  \)^60n-dM -cf-Mxdy{\-Kfs)  +  C-cf) 


where 


C  =  123 {Mx-dyf  -3U{Mx-dyf  -640 
D  =  \23(My  •  dx)4  -  384 {My  •  dx)2  -  640 


(22) 

(23) 


and  Key  -  key  /  ks ,  K/s  —kj/ks  are  the  dimensionless  effective  conductivity  and  “fluid” 
conductivity,  respectively. 

Isotropic  Laminates  Since  most  available  commercial  woven  mesh  products  are  isotropic 
structures,  i.e  dx=dy=d  and  Mx-My=M,  the  above  general  results  can  be  simplified  to  the 
following: 

cf[ 1  -  £r(«o)]  =  -3.906  y-W^nC-Md  (24) 

fi(iso)  =  4(1 -e)/d  (25) 


cf-~Ke(Jso) -  «*•***•(«:>  -l)+C-c/  .gA 


C-cf1  -  Kf 


320a  ■  Md(  1  -Kfs)  +  2  Ccf 


(26) 


Equation  24  shows  that  the  reduced  “metal  fraction”  of  isotropic  screen  laminates,  cf(]-s)  is 
solely  a  function  of  the  M/-product  of  the  mesh.  We  note  that  Md  ->  1  designates  a  “tightly” 
woven  screen.  In  actuality,  there  is  a  physical  limit  on  the  magnitude  of  Md.  For  isotropic  plain- 

weave  screens  where  the  thickness  is  limited  to  t,=2d,  M/(max)  =  1  /  ^3  =  0.577  so  that  the 
porosity  is  limited  such  that  0  <  cf(l-e)  <  0.534. 


Equation  25  shows  that  the  specific  surface  area,  |3  of  isotropic  structure  follows  the  functional 
form  of  other  porous  media,  with  SF  =  4/d.  By  contrast,  SF(spheres)  =  6/ds  for  an 
unconsolidated  bed  of  spherical  particles.  Finally,  Equation  26  shows  that  the  effective 
conductivity  Key{iso)  is  solely  a  function  of  cf,  Kjs  and  Md. 


Discussion  Figure  10  plots  the  reduced  metal  fraction  as  a  function  of  Mxdy.  The  case  for  an 
isotropic  screen  laminate  is  shown  as  a  solid  line.  This  case  always  produces  screen  laminates 
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Fig.  10  Orthogonal  weave  screen  laminate  Fig.  11  Orthogonal  weave  screen  laminate 
reduced  metal  fraction.  specific  surface  area. 

having  the  highest  reduced  metal  fraction  (lowest  porosity).  The  reduced  metal  fraction  for  two 
anisotropic  plain-weave  laminate  cases  (  M  <  1 ,  d  <  1 )  are  also  shown  in  the  figure. 

Figure  1 1  plots  (My  as  a  function  of  reduced  metal  fraction.  The  case  for  an  isotropic  screen 
laminate  is  shown  as  a  bold  solid  line.  The  shape  factor  for  this  condition  is  SF  =  4  (the  slope  of 
the  plotted  line).  This  is  true  for  all  conditions  where  d  =  1 .  The  case  M  =  0.5 ,  d  =  1  is  shown 
superimposed  on  the  isotropic  case.  All  conditions  with  M<1,  d  <1  lead  to  larger  values  of  (My 
at  fixed  Mxdy.  Therefore,  the  condition  d  =  1  always  produces  the  smallest  specific  surface  area. 

Figure  1 1  compares  (My  of  screen  laminates  with  a  bed  of  fused  spheres  and  metal  foams.  The 
theoretical  porosity  of  a  packed  bed  of  unconsolidated  spheres  depends  on  the  packing 
arrangement  [Kaviany,  1995],  It  can  range  from  e  =  0.26  for  face  centered  cubic  packing  to  e  = 
0.476  for  simple  cubic.  However,  the  achievable  porosity  of  a  packed  bed  is  difficult  to  control.  A 
typical  porosity  for  an  unconsolidated  bed  will  range  from  0.37  to  0.41.  As  a  consequence,  the 
metal  fraction  can  range  from  0.59  to  0.63  and  (M  will  then  range  from  3.54  to  3.78.  The  range  of 
the  specific  surface  area  achievable  with  uncompressed  metallic  foams  [Ashby  et  al,  2000]  is  also 
displayed.  Depending  on  the  origin  of  the  data  for  foams,  screen  laminates  provide  either  a 
significant  increase  in  specific  surface  area,  or  they  are  comparable  to  the  foams.  In  any  case, 
screen  laminates  extend  the  range  of  metal  fraction  (or  porosity)  beyond  what  appears  to  be 
achievable  with  metal  foam. 

Figure  12  plots  the  reduced  effective  conductivity  as  a  function  of  reduced  metal  fraction.  The 
case  for  an  isotropic  screen  laminate  with  Kfs  =  0  is  shown  as  a  bold  solid  line.  An  additional 

case,  isotropic  screen  laminate  with  Kfs=0.\  is  also  shown  (solid  line).  Under  these  conditions, 
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the  effective  conductivity  is  shifted  upward 
approximately  proportional  to  the  value  of 
Kfs .  In  many  applications  where  the  second 

phase  material  is  a  hydrocarbon  with  the  first 
phase  a  “good”  conductor,  Kfs<  0.005 .  Two 

additional,  anisotropic  screen  laminate  cases 
are  shown  in  the  figure.  Cases  with  M  <  1  or 
d  <  1  always  produce  effective  conductivities 
that  are  greater  than  obtained  with  isotropic 
plain-weaves. 

It  is  easy  to  find  the  maximum  possible  value  of  Key .  When  Mxdy  -» 1 ,  M  -»  Oand  d  0 ,  then, 
Eq.  (21)  gives, 

to  (max)  =  —  ^kfs  +  ^  (27). 

4  4 

The  dimensionless  effective  thermal  conductivity  of  anisotropic  screen  laminate  structures  could 
be  easily  configured  with  a  value  of  more  than  0.3,  with  the  reduced  metal  fraction  of  about  0.4. 
Comparing  with  a  fused  bed  of  spheres,  which  are  expected  to  have  Ae/spheres)  «  0.15  with 
metal  fraction,  1-e  «  0.6,  this  represents  an  approximate  two-fold  increase  in  thermal 
performance  relative  to  fused  sphere  heat  transfer  matrices.  The  value  ranges  of  effective 
conductivity  and  metal  fraction  of  the  metal  foam  [Calmidi  and  Mahajian,  1999]  are  also  shown 
in  the  figure.  By  virtue  of  their  higher  effective  thermal  conductivity,  the  screen  laminates 
structures  offer  a  two-  to  four-fold  improvement  in  performance  relative  to  high-porosity 
exchange  matrices  fabricated  from  metallic  foam. 

Model  Verification.  Transient  thermal  diffusivity  measurements  of  screen  laminate  samples 
(Appendix  A)  shows  that  the  prediction  of  the  present  model  is  within  10%  of  all  of  the 
measurements. 

Orthogonal-Weave  Heat  Transfer  and  Pressure  Drop  Correlation 

Tong  and  London  [1957]  reported  measurements  of  friction  factor  and  mesh  heat  transfer 
coefficient  for  inline  plain  weave  laminates  and  staggered  cross-rod  matrices  (no  interweaving). 
They  used  a  calorimetric  method  to  measure  the  local  heat  transfer  coefficient  of  one  heated 
filament  inside  the  array.  Since  other  filaments  upstream  of  the  “measurement  filament”  were  not 
heated,  their  correlations  are  expected  to  predict  heat  transfer  coefficients  that  are  higher  than 
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would  be  expected  if  all  wire  filaments  of  the  array  were  heated.  Miyabe  et  al.  [1982]  report  heat 
transfer  coefficient  correlations  for  plain-weave  screen  laminates  that  are  in  close  agreement  with 
those  of  Tong  and  London.  However,  they  do  not  describe  their  measurement  techniques. 
Armour  and  Cannon  [1968]  report  pressure  drop  correlations  for  plain-weave  screens,  but  not 
laminations  of  screens. 

We  postulate  that  the  pressure  drop  across  the  mesh,  AP,  is  functionally  related  to  fluid  flow 
properties  and  divergence  angle  d>  (Fig.  2)  as  follows: 

AP 

—  =  (28) 
where  G  =  —  is  the  internal  mass  velocity.  Then  dimensional  analysis  gives 


/  =  ./h(Re,£,<D) 


(29) 


where  f  = 


2  pAP 
G  2  fit 


is  the  friction  factor,  and  Re  = 


GDh 


is  the  mesh  Reynolds  number,  with 


Dh  = 


4s 

1 


sd 

!-£• 


the  mesh  hydraulic  diameter. 


In  a  similar  way,  we  postulate  that  the  mesh  heat  transfer  coefficient,  h  is  functionally  related  to 
fluid  and  flow  properties  as  follows: 

h  =  f{p,G,n,d,p,t,e,  Pr,0)  (30) 

Dimensional  analysis  gives 

St  =  fh(Re,  Pr ,d  /  t,e,< I>)  (31) 

where  Pr  is  the  Prandtl  number  of  the  coolant.  In  the  following,  we  describe  experiments  to 
determine  the  specific  form  of  eqs.  (29)  and  (31). 

Heat  Transfer  and  Pressure-Drop  Correlations,  Orthogonal  Flow  Experiments  are  performed 
to  measure  the  pressure  drop  and  porous  wall  effective  conductance  for  isotropic  screen  laminates 
with  the  approach  flow  normal  to  the  screen  lamination  (O  =  90°).  Then,  the  mesh  Stanton 
number  is  determined  from  eq.  (8).  The  friction  factor  and  Stanton  number  correlations  are  found 
in  the  case  of  inline  and  staggered  configuration  respectively.  The  test  faciloity  and  experimental 
procedures  are  described  in  Appendix  B. 
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Figure  13  summarizes  pressure  drop 
measurements  for  inline  and  staggered  screen- 
laminates.  The  staggered  screen-laminate 
friction  factor  is  much  higher  than  the  friction 
factor  of  the  inline  configuration.  The  present 
data  for  the  inline  correlation  is  correlated 
with  the  following  expression 


/  = 


0.39  + 


24.5 

Re 


v  0.4 


1  -6j 


(32) 


Fig.  13  Friction  factor  of  inline  and  staggered 


.  .  ,  .  .  .  star.kftH  sr.rpip.n  laminates 

The  expression  contains  an  inertial  loss  term 


and  a  viscous  term.  We  note  that  - 

\-e 


These  results  are  compared  with  Tong  and 


London’s  data.  The  present  screen  laminates  are  bonded  whereas  those  of  Tong  and  London  are 
not;  bonding  fillets  at  wire  intersections  probably  gives  rise  to  higher  inertial  losses  at  high 
Reynolds  numbers.  Eq.  (32)  reproduces  the  data  that  generated  it  with  a  standard  error  of  ±20%. 


Figure  13  also  shows  the  friction  factor  vs.  Reynolds  number  correlation  for  staggered  screen 
laminate  stacking.  In  this  case,  the  data  is  correlated  as  follows 


55.3 

Re 


(33) 


The  data  show  that  staggered  stacking  results  in  an  approximate  two-fold  higher  pressure  loss 
than  inline  stacking  at  the  same  flow  rate.  Eq.  (33)  reproduces  the  data  that  generated  it  with  a 
standard  error  of +15%. 


Heat  transfer  experiments  are  performed  with  chilled  water  (Pr  «  9.5).  We  are  correlating  data  in 
terms  of  Colburn  j-factor  so  that  a  direct  comparison  can  be  made  with  the  data  reported  by  Tong 
and  London,  who  experimented  with  air  (Pr  =  0.7).  We  have  chosen  to  correlate  the  data  with  a 
power-law.  A  regression  analysis  yielded  the  following  results. 


j  =  St  -  Pr3  =  5.86Re~ 


0.68 

1 

to 

f£ 

K\-e)  1 

(in-line) 


(34) 


j  =  St-  Pr3  =2.63  Re 


-0.65  [  G 


.-0.2 


l-e y 


f£ 


(staggered)  (35) 
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Figure  14  summarizes  heat  transfer  data  (both 
inline  and  staggered  configurations)  in  terms 
of  modified  j-factor  as  a  function  of  Reynolds 
number.  The  modified  j-factor  correlation  for 
the  inline  correlation  is  plotted  as  continuous 
line  and  compared  with  Tong  and  London’s 
correlation  for  air.  Tong  and  London’s  data 
are  adjusted  to  conform  to  the  definition  in  the 
present  case  and  then  plotted  as  a  dashed  line. 

Fig.  14  Modified  Colburn  j-factor  of  inline  and 

Comparing  the  present  correlation  with  the 

Tong  and  London’s  data,  the  magnitude  of  the  present  data  is  smaller  than  the  correlations  of 
Tong  and  London  at  the  higher  Reynolds  number.  Otherwise,  the  results  are  in  close  agreement. 
This  implies  that  the  present  correlation  is  applicable  to  coolants  with  0.7<Pr<9.5.  The  figure 
also  shows  that  the  staggered  configuration  produces  a  lower  modified  j-factor  than  does  the 
inline  configuration.  This  is  a  somewhat  surprising  result.  Eqs.  (34,  35)  reproduce  the  data  that 
generated  them  with  a  standard  error  of  ±18%. 


0.57  <£-<0.84 

"  ♦  '  Present  (inline) 

_  A Pres  ert(stagg  ere  d) 

-  -  -  Tong  and  London[1957], 
Screen  rr=1 

*  * 

% 

> 

Heat  Transfer  and  Pressure-Drop  Correlations,  Non-Orthogonal  Flow  Both  pressure  drop  and 

heat  transfer  are  measured  for  a  single  layer  of  isotropic  screen  with  air  as  described  in  Appendix 

B.  Test  articles  with  different  divergence  angles  (6°  <2<p  <  180°)  and  laminate  heights  (16.1mm 

<  H  <  28.4mm)  are  investigated.  All  of  our  test  articles  are  made  of  solder-bonded  copper  screen 

with  d  =  0.48  mm  (0.019”),  M  =  7.9  cm'1  (20in_1)  with  e  =  0.677,  ft  =  2.68x1 03/w_I,  and 

kev  _  55  2  wa“  as  estimated  by  Eq.  (24  -  26) 
m-K 


In  Fig.  15,  we  present  just  several  sets  of  pressure  drop  data  (<p  =3°,  7.1°,  20°,  34°,  90°)  in  terms 


of  modified  friction  factor,  F  =  f  ■ 
expression: 


The  data  is  best  correlated  with  the  following 


( 


F  = 


0.39  + 


24.5 

Re  , 


+ 


0.17 


26.3 


tan(^)  Re  tan(^>) 


(36) 


The  expression  contains  both  inertial  loss  terms  and  viscous  terms.  The  first  two  terms  are  for  <p 
=90°  (Eq.  32),  and  the  second  two  terms  reflect  the  effect  of  variations  in  screen-fin  divergence 
angle.  Eq.  (36)  reproduces  the  data  that  generated  it  with  a  standard  error  less  than  ±10%. 
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Fig.  15  Modified  friction  factor  non-orthogonal  Fig.  16  F  factor  versus  divergence  angle 
flow. 


Fig.  17.  J  Factor  of  Screen  Laminate  Heat  Sink  Fig.  18.  Non-orthogonal  flow  modified  J  Factor 

versus  cp,  Eq.  (37) 

Figure  16  plots  the  modified  friction  factor,  Eq.(36)  as  a  function  of  <p  for  three  Reynolds 
numbers.  As  shown  in  the  figure,  F  decreases  very  sharply  when  0  <  0  <  20°  and  the  curve 
becomes  flat  when  the  0  is  greater  than  about  40°. 


Figure  17  summarizes  heat  transfer  data  in  terms  of  modified  Colburn  j-factor, 


J  -  j ' 


\-s 


—  as  a  function  of  Reynolds  number.  Data  for  <p  =  3°,  6.2°,  9.5°  and  cp  >  20° 


*  A d, 

are  plotted.  The  figure  shows  that  J  is  sensitive  to  <p  when  it’s  less  than  about  14°;  and,  variations 
in  0  have  no  significant  effect  on  J  when  0  is  greater  than  14°.  We  have  chosen  to  correlate  the 
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data  with  a  power-law.  A  regression  analysis  yielded  the  following  results.  Eq.  (37)  reproduces 
the  data  that  generated  it  with  a  standard  error  of  ±7  %. 


J  = 


A(<p)  Re" 72  0°<<p<  8° 

4.6[A(<p  - 4.3°)]°9  -  2.12[A(<p  -  4.3°)} 2  Re'072  8°  <  q>  <  14° 
2.36  Re-0 72  14 °  <<p 


(37) 


where  A(<p)  =  104s“(?,)-03 . 

In  Figure  18,  we  plot  J  (Eq.37)  as  a  function  of  (p  for  three  flow  rates.  The  figure  shows  the  rapid 
increase  in  J  with  increase  in  divergence  angle  up  to  q>  »  14°,  with  J  constant  for  <p  >  14°. 

Multiple-Row  Heat  Transfer  &  Pressure  Drop  A  single  serpentine  row  of  screen  laminate 
cannot  provide  enough  surface  area.  In  order  to  gain  more  surface  area  we  need  to  pattern  several 
rows  of  screen  laminate  in  a  series  arrangement.  Based  on  the  test  data  for  multiple  rows,  we 
develop  friction  factor  and  Colburn  j-factor  corrections  that  account  for  the  decrease  in  row-by- 
row  performance  of  a  multiple  row  system.  Here  we  define  the  overall  pressure  drop  and 
conductance  as  AP„  and  U(tW)  as 


^*nrow  ^  row  1  —  —  6 

(38) 

(39) 

From  the  above  equations,  it  is  evident  that, 
with  increases  in  the  number  of  screen-fin 
rows,  the  overall  conductance  increases  more 
rapidly  than  does  the  overall  pressure  drop. 
So,  a  larger  number  of  rows  is  favored  with  a 
design  that  seeks  to  optimize  performance  of  a 
heat  sink  having  a  pressure  drop  constraint. 

Equations  (36)-(39)  can  be  used  to  design 
screen-laminate  heat  sinks.  In  order  to  test  the 
accuracy  of  the  design  algorithm,  we 
constructed  a  prototype,  measured  its  pressure 
drop  and  heat  transfer  performance,  and 
compared  these  results  with  performance 
predictions  embodied  in  above  equations. 


Fig.  19  Prototype  Heat  Sink 
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Table  1  Comparisons  Between  Test  Data  and  Model,  20  =  22°. 


Conductance  |watt/°C] 

4  rows 

5  rows 

6  rows 

Model 

Test  Data 

Model 

Test  Data 

Model 

Test  Data 

24.907 

1.15 

1.01 

1.22 

1.105 

1.205 

1.17 

AP[Pa] 

62.267 

1.37 

1.255 

1.51 

1.465 

1.555 

1.535 

99.628 

1.47 

1.37 

1.66 

1.645 

1.695 

Standard  Error 

10.5% 

6.3% 

2.5% 

The  prototype  is  a  simple  finned  heat  sink,  63.5  mm  deep,  45.7  mm  wide  with  23.4  mm  high 
screen-fins  {d  =  0.483mm,  M  =  787.402  m1),  as  shown  as  Figure  19.  It  consists  6  rows  of 
serpentine-pattern  screen-fins  with  divergence  angle,  2<p=  22°.  We  also  studied  how  the  row 
number  affects  performance  by  successively  cutting  the  last  row  from  the  piece  after  each 
experiment  so  that  we  measured  the  performance  with  4  rows,  5  rows  and  6  rows.  In  Table  1  we 
present  comparisons  between  the  test  data  and  prototype.  The  table  shows  that  our  model  agrees 
well  with  our  test  data.  The  standard  error  is  less  than  11%. 

Non-orthogonal  Weave  Porosity,  Specific  Surface  Area,  Effective  Thermal  Conductivity 
Figure  20  summarizes  the  geometric  parameters  of  the 

physical  model  of  a  diamond-shape  plain-weave  screen. 

The  figure  shows  a  segment  of  a  plain-weave  that 

comprise  a  unit  cell  of  the  material,  together  with  a 

section  view  that  illustrates  the  weave  pattern  of  the 

wire  filaments.  Wire  filament  axes  are  directed  at  #, 

and  #2  with  respect  to  the  presumed  heat  flux  direction 
(q-direction),  so  the  weave  angle  is#,  +  #2  -26. 

Serpentine  wire  filaments  have  diameter  di  and  d2  and 
corresponding  mesh  numbers  ML  and  M2. If  we  consider 
symmetric  weaves  only,  then  the  mesh  number 
M=Mi=M2,  filament  diameter  d=di=d2  and  the 
directions  of  wire  filament  axes  relative  to  heat  flux 
direction#,  =#2  =  #. 

Fig.  20  Diamond-weave  unit  cell. 

Fig.21  shows  a  typical  diamond-shape  solder-bonded 

copper  plain-weave  screen  laminate.  The  wire  filaments  have  an  effective  diameter  of  0.48mm, 


4  ! 


i 

A-A  Section 


and  the  mesh  number  is  6.3cm'1  (16  inch"1). 
The  sample  has  s  =  0.724.  (After  the 
bonding,  the  effective  thickness  of  solder  layer 
over  the  screen  filament  ranges  from  0.006mm 
to  0.013mm.) 

Consideration  of  Fig.  20,  shows  the  ratio  of 
the  reduced  metal  fraction  of  a  diamond- 
weave  to  an  equivalent  orthogonal  weave  can 
be  written  as  a  function  of  the  weave  angle 
(2 6)  only;  and  the  same  result  is  obtained 
with  respect  to  the  cff}- product. 

tfnjl  '(1  ~g<r/2)  _  tfxjl  '  Ar/2  _  gjjjO Q\ 

cfie  '  (1 —  eie  )  cfx?  ’  Pie 


Fig.21  Bonded  Diamond-shaped  Isotropic 
Screen  Laminates 


(40) 


where  the  subscript  “n/2”  stands  for  orthogonal,  the  subscript  “ 29 ”  stands  for  diamond-shape 
with  arbitrary  weave  angle  29 .  Equation  (40)  shows  that  diamond-weaves  possess  greater  metal 
fraction  and  specific  surface  area  than  orthogonal  weaves. 


The  unit  cell  shown  in  Fig.20  is  transformed  to 
the  layered,  rectangular  cross  section  segment 
shown  in  the  Fig.22.  Each  rectangular  wire 
filament  has  thickness  ndlAg  and  width  gd  so 
that  the  cross  section  area  of  each  filament 
is  mi2 /4.  The  geometric  index,  g  is  a  measure 
of  the  contact  between  wire  filaments  at 
intersections.  We  further  require  that  the 
volume  of  wire  filaments  be  preserved  across 
the  geometric  transformation.  Then 
g  =  n/(4-cf). 

If  we  neglect  the  contribution  of  the  second 
phase  (fluid)  to  the  effective  conductivity  of 
the  weave,  then  the  in-plane  effective  thermal 


Fig.22  Diamond-weave  transformed  unit  cell 
and  thermal  circuit 


conductivity  in  the  q-direction  may  be  determined  by  considering  the  thermal  circuit  for 
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conduction  across  the  transformed  unit  cell,  which  is  also  shown  in  Fig.  22.  It  shows  four 
resistance  elements  for  conduction  along  the  axes  of  the  wire  filaments  that  make  up  the  heat  flow 
path  across  the  transformed  unit  cell.  For  a  diamond-shaped  symmetric  screen  laminate,  the 
thermal  resistance  of  each  filament  segment  is  the  same.  It  is  easy  to  show  the  total  thermal 
resistance  of  the  unit  cell  in  the  q-direction  is 


total 


-6.25-10  ~3-C 
x -ks-d2-M 


(41) 


where  C  is  given  in  Eq.  (22)  and  ks  is  the  thermal  conductivity  of  the  filament  material.  The  q- 
direction  in-plane  effective  thermal  conductance  across  the  transformed  unit  cell  is. 


jy  cf-Widx+dy ) 

K total  ~  e  '  £ 


(42) 


where  L  is  the  length  of  unit  cell  along  the  q-direction,  W  is  the  width  of  the  unit  cell 
perpendicular  to  the  q-direction,  and  ke  is  the  in-plane  effective  thermal  conductivity  along  q- 
direction.  Rearranging  Eq.  42  gives 


cf-Ke  = 


-Wn-Md 
tan  (0)-C 


(43) 


where  Ke  -  keq  /  ks ,  is  the  dimensionless  effective  thermal  conductivity  in  the  q-direction.  With 


the  weave  angle  equal  to  nil,  Eq.43  gives  the 
same  result  for  the  effective  thermal 
conductivity  of  the  orthogonal  weaves  (Eq.  26). 

We  find  that  for  orthogonal  isotropic  screen 
laminates,  that  changing  the  orientation  with 
respect  to  the  heat  flux  vectors  (i.e.  in-plane 
rotation  of  the  laminates),  will  not  change  the  in- 
plane  effective  conductivity  of  the  structure.  So 
with  the  same  filament  diameter  and  mesh 
number,  the  two  configurations  of  orthogonal 
isotropic  screen  laminates  in  the  Fig.23  have  the 
same  effective  in-plane  effective  thermal 
conductivity  in  the  q-direction. 

The  ratio  of  reduced  conductivity  of  diamond 


I 


Configuration  2 


Heat  flux  along  Heat  flux  not  along 

Fig.23  Effect  of  in-plane  rotation  on  Ke<,  of 
orthogonal  plain-weave  laminates 
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weaves  to  orthogonal  weaves  with  the  same 
mesh  number  and  filament  diameter  can  be 
written  as: 

cfie  •  Ke2o  _  sin(2fl) 

Cfic!2  '  ^e,T  /  2  1 _  COs(20) 

(44) 

The  ratio  of  conductivity  of  diamond  weaves 
to  orthogonal  weaves  is  plotted  in  Fig.24  as  a 
function  of  20.  We  note  that  due  to  the 
physical  limitation  of  the  geometry,  there  are 
upper  and  lower  limiting  values  of  the 
arbitrary  weave  angle.  For  the  symmetric 
diamond-shape  weave: 

2arcsin(M//2)<20<2arccos(Mf/2)  (45) 

Benchmark  measurement  of  the  thermal  diffiisivity  of  diamond-weave  samples  (Appendix  A) 
shows  that  the  prediction  of  the  present  model  is  within  8%  of  all  of  the  measurements. 

Exchange  Matrix  Performance  Evaluation 

Screen-laminates  offer  considerable  design  flexibility.  Adjustment  of  wire  diameters  and  wire 
pitch  allows  for  control  of  the  structure’s  porosity,  heat  transfer  surface  area-to-volume  ratio,  and 
effective  thermal  conductivity.  However,  the  friction  factors  and  Stanton  numbers  of  the  heat 
exchange  matrix  are  comparable  to  those  of  other  heat  transfer  surfaces.  The  question  that  must 
be  addressed  is:  under  what  conditions  does  the  screen-laminate  technology  offer  superior 
performance. 

In  the  following,  we  describe  a  fixed  outer  geometry  comparison  [Webb,  1994]  of  the 
performance  of  a  screen-laminate  exchange  matrix  with  an  exchange  matrix  consisting  of  an 
unconsolidated  bed  of  spherical  particles  having  the  same  mass  (porosity)  and  face  area.  The 
particle  bed  is  selected  since  both  exchange  matrices  can  then  be  treated  as  porous  media,  so  that 
the  thin-fin  model  can  be  used  to  predict  overall  thermal  performance.  In  this  way,  performance 
model  selection  will  not  influence  the  ranking  of  performance. 


Fig.24  The  ratio  of  conductivity  of  diamond- 
weave  to  orthogonal  configuration. 
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Table  2  Screen  laminate  and  sphere-bed  characteristics 


Characteristic 

Screen-Laminate 

Spheres 

External  Characteristics 

£=0.4,  \0<  t/d  <100, 10^  H/d <100 

Copper  solid  phase  (ks=  400  W/mK) 

Coolant  is  air  (a)y  300K  (Pr  =  0.69) 

Internal  Geometric 

d  =  0.5mm  ,  Md= 0.556, 
c/=0.85 

ds  =  0.5  mm 

Effective  Conductivity 

88  W/mK,  Xu  and  Wirtz15 

34  W/mK,  Hadley4 

Specific  Surface  Area 

4800m*1,  Eq.  (25) 

7200  m*1,  Kaviany5 

Pressure  Drop 

Eq.  (32) 

Dullien3 

Matrix  Heat  Transfer 
Coefficient 

Eq.  (34) 

Wakao  and  Kaguei12 

The  thermal  performance  of  a  plain-weave 
screen  laminate  having  porosity  (s  =  0.4)  is 
compared  with  a  spherical  bed  system  of  the 
same  porosity.  This  can  be  achieved  with  a 
screen-laminate  system  having  Md=  0.556 
with  c^=0.85.  This  will  result  in  the  two 
systems  having  the  same  mass  and  external 
dimensions. 

The  details  of  the  two  systems  are 
summarized  in  Table  2.  Computations  are 
performed  with  air  as  coolant.  Other,  higher 
Prandtl  number  fluids  produce  similar  results. 

External  dimensions,  thickness  and  height  are 
varied  over  a  wide  range,  10  <  t/d  <  100,  10  < 

H/d  <  1 00,  while  the  screen  laminate  wire  diameter  and  the  sphere  diameter  are  held  equal  to  each 
other,  d  =  0.5  mm.  Calculations  for  other  diameters  give  essentially  the  same  results.  The  table 
shows  that  the  screen  laminate  system  possesses  a  higher  effective  thermal  conductivity  by  a 
factor  of  about  2.6,  whereas  the  sphere  bed  has  a  higher  specific  surface  area  (by  a  3:2  ratio). 


Fig.  25  NTU  and  pressure-drop  performance 
ratios. 


We  define  the  ntu-ratio  and  pressure  drop  performance  ratios  as  follows: 


ntu(screens)  _  A P(screens) 

ntu(spheres)  p~srad  AP(spheres) 


(46a,b) 
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Fig.  26  Screen  laminate-to-sphere  bed  Fig.  27  Screen  laminate-to-sphere  bed  thermal 
thermal  performance  ratio,  Ry  at  equal  performance  ratio,  Ru  at  equal  pressure  drop 
superficial  mass  velocity 


Figure  25  plots  these  ratios  as  a  function  of  Reynolds  number  based  on  filament  size, 
Rerf  =G,dl  p  for  exchange  matrices  with  t/d  =10  (the  base  line  configuration).  Red  reflects 

variation  in  superficial  mass  velocity  incident  on  the  exchange  matrix.  The  figure  shows  that  at 
equal  superficial  mass  velocity,  the  sphere  bed  number  of  transfer  units  is  four  to  six  times  greater 
than  that  of  the  screen  laminate  system.  This  is  due  to  the  higher  heat  transfer  coefficient  (h)  and 
specific  surface  area  (P)  exhibited  by  the  sphere  bed.  On  the  other  hand,  the  pressure  drop  across 
the  sphere  bed  is  seven  to  ten  times  greater  than  that  across  the  screen  laminate. 


The  results  of  Fig  25  imply  that  sphere  beds  should  outperform  screen  laminate  systems  at  equal 
coolant  flow  rates.  However,  the  above  analysis  does  not  reflect  the  2.6-fold  higher  effective 
thermal  conductivity  of  the  screen  laminate  system.  In  order  to  assess  this  component  of  the 
problem,  we  define  the  “heat  duty  ratio” 


U {screens) 
U  {spheres) 


(47) 


where  U,  defined  by  Eq.  (1),  is  calculated  from  Eq.  (8).  Figure  26  plots  Ry  as  a  function  Rea. 
The  base-line  configuration  {t/d  =  10,  H/d  =  10)  shows  the  sphere  bed  outperforming  the  screen- 
laminate  system  by  up  to  a  factor  of  two  at  high  coolant  flow  rates.  However,  an  increase  in 
exchange  matrix  height  {H/d  =  100)  gives  rise  to  a  profound  change  in  results.  In  this  case,  the 
screen-laminate  system  out-performs  the  sphere  bed  by  up  to  40%.  This  is  a  reflection  of  the 
higher  effective  thermal  conductivity  of  the  screen-laminate  system.  On  the  other  hand,  an 
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increase  in  exchange  matrix  thickness  (relative  to  the  base  line  configuration)  gives  rise  to  a  slight 
decrease  in  performance.  It  should  be  noted  that  while  Eqs.  (32,  34)  include  the  effect  of 
laminate  thickness  on  mesh  heat  transfer  coefficient,  the  Wakao  and  Kaguei  Nusselt-number 
correlation  for  sphere  beds  does  not. 

The  situation  is  further  improved  (relative  to  sphere  bed  thermal  performance)  if  the  comparison 
is  done  on  the  basis  of  equal  pressure  drop.  Figure  27  plots  Ru  versus  the  dimensionless  pressure 
drop,  fRerf2.  In  this  case  the  base-line  configuration  (t/d  =  10,  H/d  =  10)  outperforms  the  sphere 
bed  by  as  much  as  a  factor  of  7.  However,  it  does  so  at  he  expense  of  vastly  increased  coolant 
flow  rate,  as  reflected  by  a  plot  of  the  mass  velocity  ratio, 

_  Gi(screens) 

G‘  Gi(spheres ) 

where  the  screen  laminate  flow  rate  is  seen  to  be  as  high  as  7-times  that  for  the  sphere  bed.  An 
increase  in  exchange  matrix  height  gives  rise  to  a  reduction  in  performance  advantage  at  low 
pressure  drop  and  a  slight  increase  in  performance  advantage  at  high  pressure  drop. 

A  P4-Format  Heat  Sink  Based  on  Screen  Laminate  Technology 

The  previous  analysis  indicates  that  air  cooled 

heat  sinks  based  on  screen-laminate 
technology  are  promising  altumatives  to  dense 
folded-fin  heat  sinks  that  are  currently  used  in 
commercial  microprocessor  cooling 
applications.  In  the  following,  we  describe  a 
novel  high-performance  heat  sink  based  on 
screen-fin  technology.  The  objective  is  to 
design  a  P4-format  air  cooled  heat  sink, 
having  a  3”  x  2.5”  footprint  and  2”  height  that 
will  have  an  overall  thermal  resistance 
(excluding  the  base  plate  spreading  resistance)  Fig-  28  Schematic  of  screen-fin  heat  sink 
of  0.2  °C/watt  at  0.25”  H20  pressure  drop  of  the  coolant.  This  is  the  unofficial  industry 
benchmark  for  a  longitudinal  fin  heat  sink. 

Screen  Laminate  Heat  Sink  Design  Concepts  Figure  28  shows  top  and  side  views  of  a  screen- 
fin  heat  sink.  The  heat  sink  consists  of  a  base  plate  (presumed  to  be  in  thermal  contact  with  the 
electronics)  and  a  fin  structure.  The  fin  structure  consists  of  a  screen-fin  oriented  perpendicular 
to  the  heat  sink  base,  and  laid  out  in  a  serpentine  arrangement  so  that  the  divergence/convergence 


GoJ. 


f,o 
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angle  between  segments  of  the  screen  is  2<p.  The  external  envelope  is  Lo  wide,  W0  deep,  and  Htot 
high.  H=Htoi  -  Hb  is  the  fin  height.  A P  is  the  pressure  drop  across  the  serpentine  screen-fin  heat 
sink.  Go,  Tfo  denotes  the  inlet  air  mass  velocity  and  temperature  respectively. 

As  shown  in  the  figure,  the  screen-fin  is  laid  out  in  a  serpentine  arrangement  so  that  the  ambient 
air  with  uniform  temperature  T/i0  and  flow  rate  Go  flows  down  a  converging  section,  through  the 
screen  and  then  out  through  a  diverging  section.  This  arrangement  provides  for  an  increase  in 
heat  transfer  surface  area,  and  it  reduces  the  superficial  mass  velocity  incident  on  the  screen 
laminate,  thereby  effectively  controlling  coolant  pressure  drop. 


Consider  the  one-row  heat  sink  segment 
shown  in  Fig  29.  The  heat  sink  has  depth  Wo 
=  63.5mm,  width  Lo  =45. 72mm,  and  fin 
height,  H  =  38.1mm.  At  fixed  volume  flow 
rate,  Q  =  G0(LoH)l p,  variations  in  6  have  two 
effects  on  the  pressure  drop  and  overall 
conductance: 

1.  Increasing  <p  will  result  in  a  decrease 
in  heat  transfer  surface  area  (As)  and 
friction  factor.  However  the  mesh 
Stanton  number  will  increase  as  per 
Eq.  (37)  and  mesh  Stanton  number 
will  vary  with  <p. 

2.  On  the  other  hand,  the  heat  transfer 
surface  area  will  increase  with 
decreasing  <p  since  more  screen-fin 
segments  can  be  placed  on  the  base 
plate.  However,  at  fixed  flow  rate,  an 
increase  in  surface  area  will  result  in  a 
decrease  in  G„  resulting  in  a  decrease 
in  Stanton  number. 

Figure  30  shows  the  overall  effect  of  variation 
in  (p  on  overall  conductance  of  the  segment 
when  Q  =  9.439  liter/sec.  The  figure  shows 


Fig  29  One-row  heat  sink  segment 


Divergence  Half-Angle  [degree] 
Fig.  30  Effect  of  Angle  of  incidence 
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that  8°  <  <p  <13°  produces  the  highest  unit 
surface  conductance  (depending  on  the 
number  of  screen-fin  rows),  with  lower 
conductance  indicated  at  higher  and  lower 
values  of  <p.  This  analysis  also  shows  that  the 
overall  pressure  drop  is  minimum  when  11°  < 

<p  <  20°.  Approximately  the  same  result 
obtains  for  other  flow  rates. 

Figure  31  shows  that  the  conductance  of  the 
screen-fin  heat  sink  segment  changes  with  row 
number  with  other  conditions  fixed.  At  small 
nmw,  the  surface  area  increases  with  the  row 
number,  so  does  the  performance  of  the  heat  sink.  But  with  further  increase  in  nmw  the  flow  rate 
reduces,  leading  to  a  decrease  in  overall  performance.  For  given  Md-product  and  applied  AP 
there  will  be  an  optimum  value  of  nrow.  The  figure  shows  that  for  the  current  overall 
configuration  with  6-11°  AP  —  62  Pa  (0.25  inch  H2O),  an  optimum  exists  at  nrow  =  5,  Md  = 
0.528.  Other  optimums  obtain  for  other  configurations. 

Table  3  Screen-fin  heat  sink  characteristics.  L0  =  63.5mm,  We  =  76.2mm,  H  =  38.1mm 


Parameter 

Screen-Fin  Heat  Sink 

Material 

Copper  Plain-weave 

Fin  spec 

d  =  0.56mm  (0.022”),  M  =  9.45  cm’1  (24m  \ 

e 

14° 

Thermal  Conductivity 

72  watt/(m*K) 

Brow 

5 

H.T.  Surface  Area 

0.154  m2 

Weight  of  Fin  Structure 

193  gram 

Coolant  flow  (2}  AP  =  62  Pa 

8.0  liter/sec 

Optimum  Design  of  Screen  Laminate  Heat  Sink  Based  on  the  above  discussion,  an  optimum 
design  does  exist  for  a  given  size  (L„,  Wa  H)  and  applied  pressure  drop,  AP„.  If  we  consider  a 
heat  sink  with  a  76.2mm  x  63.5mm  (3”  x  2.5”  footprint)  having  38.1mm  (1.5”)  tall  screen-fins 
and  operating  at  AP  =  62  Pa  (0.25”  H20),  then  systematic  variation  of  design  parameters  results  in 
the  optimal  design  copper  heat  sink  with  properties  summarized  in  of  Table  3.  The  fin  format  of 
this  prototype  is  as  shown  in  Fig.  19.  This  heat  sink  will  have  an  overall  conductance  of  4.3 
watt/°C  (thermal  resistance  =  0.23  °C/watt)  when  operated  at  APn  =  62.27  Pa.  This  is  within  14% 
of  the  target  benchmark,  well  within  the  error  tolerance  that  we  would  place  on  the  analysis. 


Fig.  31  Effect  of  Number  of  Rows,  AP  =  62  Pa. 
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Assuming  a  0.25”  thick  base  plate  (typical  for 
this  application)  the  prototype  is  0.25”  shorter 
than  the  target.  A  conventional  heat  sink 
would  of  necessity  be  of  folded  fin  design, 
with  a  very  fine  fin  pitch.  On  the  other  hand, 
the  screen-fin  prototype  is  quite  simple  in 
construction.  The  mass  of  the  fin  structure 
(193  gm)  is  about  the  same  as  would  be  of  an 
aluminum  folded  fin  design.  Furthermore, 
since  the  present  prototype  is  all  copper,  the 
usual  issues  that  accompany  bonding 
aluminum  fins  to  a  copper  base  plate  would 
performance  of  the  screen-fin  prototype. 


Fig.  32  Optimum  design  of  screen- 
fin  heat  sink. 


not  be  present.  Figure  32  plots  the  expected 
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Three-Dimensional  Orthogonal  Weaves 

The  following  material  demonstrates  the  fabrication  methodology  of  a  3-D,  aluminum  wire 
filament,  bonded  mesh  deployed  as  a  heat  exchange  surface.  Models  of  the  porosity,  specific 
surface  area  and  effective  thermal  conductivity  of  the  mesh  is  developed.  Mesh  Stanton  number 
and  friction  factor  correlations  for  a  coolant  with  Prandtl  number  equal  to  9.5  (chilled  water)  are 
reported.  Measurements  are  reported  for  fabricated  test  samples  of  varying  thickness.  A  heat 
exchanger  performance  evaluation  comparing  the  3D  woven  mesh  technology  to  screen  laminate 
technology  is  described. 

We  have  found  that  the  weaving/wire  bonding  process  must  be  carefully  controlled  to  insure  that 
target  porosity,  specific  surface  area  and  effective  thermal  conductivity  are  achieved.  Effective 
thermal  conductivities  are  found  to  be  at  least  two-times  larger  than  achieved  in  other  comparable 
porous  media  configurations.  Mesh  friction  factor  and  Stanton  number  are  comparable  to  those 
achieved  with  other  exchanger  surface  technologies.  An  exchanger  performance  evaluation 
shows  that  exchangers  having  superior  performance  can  be  configured. 

Stacked  Weave  Geometry 

Figure  33  shows  a  three-dimensional 

orthogonal  stacked-weave  that  consists  of 
three  separate  wire  filaments  of  diameters  dx, 
dy  and  dz,  having  axes  aligned  with  the  x,  y 
and  z-axes,  respectively.  Coolant  flow  is 
presumed  to  be  primarily  in  the  x-z  plane.  The 
y-wire  filament  diameter,  dy  is  larger  than  dx 
and  dz  so  that  the  effective  thermal 
conductivity  in  the  y-direction,  key  is  larger 
than  kex  or  kez.  In  this  way,  heat  is  transferred 
to  the  fluid  by  conduction  primarily  along  the 
y-filaments;  then  to  the  x-  and  z-filaments,  which  act  as  fins.  The  wire  filaments  are  bonded  at 
intersections  to  facilitate  conduction.  We  designate  the  three-dimensional  orthogonal  weave 
shown  in  Fig.  1  a  “stacked  weave”  since  there  is  no  interweaving  of  wire  filaments  in  any  of  the 
three  principal  planes.  This  approach  allows  for  a  very  dense  structure. 

Mesh  Fabrication 

Our  design  objective  is  to  fabricate  three  3D  stacked  weave  heat  exchange  matrixes  that  are 
approximately  20mm  high  x  100mm  wide  with  t  =  6.35mm,  9.53mm  and  12.7mm,  respectively. 


Fig.  33  Three-dimensional  orthogonal  “stacked” 
weave. 
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The  mesh  is  to  have  dy  =  0.76  mm  (0.030”), 
dx  =  dz  =  0.38mm  (0.015”).  Analysis  in  the 
following  section  shows  that  the  target  metal 
fraction  is  (1-e)  =  0.61,  and  f3  =  4581m'1. 

The  target  effective  thermal  conductivity  in 
the  y-direction  is  84  W/mK. 

The  fabrication  process  involves  weaving  a 

wire  mesh  “rope”  of  specified  wire  filament 

diameters  and  mesh  numbers  in  the  three 

coordinate  directions,  braze-bonding  the  _  ,  ,  ,  _ . 

Fig.  34  Three-D  woven/bonded  mesh  fabrication 

wire  filaments  at  their  intersection  points,  sequence, 
and  cutting  (via  the  wire-EDM  process)  and 

braze-bonding  mesh  segments  together  to  form  the  heat  exchange  matrix.  Figure  34  summarizes 
the  process.  Aluminum  wire  mesh  samples  (Alloy  1100)  were  woven  on  a  commercial  loom 
[T.E.A.M.  Inc,  slatersville  RT|.  The  large  diameter  wire  filaments  (dy)  formed  the  warp  (weave) 
direction,  and  dx-wire  and  dz-wire  filaments  made  up  the  shut  and  fill  directions.  The  loom 
restricted  the  3D  weave  to  a  product  that  was  approximately  35mm  wide  x  15mm  thick. 

The  woven  “rope”  was  dip-brazed  to  bond  wire  filaments  at  their  intersection  points.1  The 
process  involves  impregnating  the  woven  structure  with  a  slurry  of  al/si  eutectic  alloy  in  flux, 
heating  the  article  to  the  flux  activation  temperature,  and  then  immersing  the  article  in  a  salt  bath 
at  590°C. 

Since  the  dy-wire  filaments  are  in  the  warp-direction,  and  the  height  (y-direction)  of  the  mesh  test 
articles  is  20mm,  twenty  mm  long  segments  of  the  woven/bonded  rope  had  to  be  cut,  rotated  90°, 
and  then  bonded  together  to  form  the  test  article.  This  was  accomplished  via  wire-EDM. 

Table  4  summarizes  results  of  the  fabrication  process.  The  table  lists  target  and  achieved  mesh 
physical  and  thermal  attributes.  Achieved  physical  dimensions  were  determined  via  microscopic 

Table  4  Three-dimensional  mesh  characteristics. 


Parameter 

Target 

Achieved 

A% 

dx,  dz,  dy  [in] 

0.015, 0.030 

0.0153,  0.031 

+1.7 

Mx,  Mz  [1/in] 

22.2 

20.7 

-2.3 

My  [1  /in] 

33.3 

31.5 

-5.5 

e,  Porosity 

+26 

nm 

4580(eq  5) 

-16.2 

84  (eq.  13) 

75.2 

-10.5 
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inspection  of  sectioned  mesh  samples,  and  the 
porosity  was  determined  via  gravimetric 
measurements.  Details  regarding  these 
measurement  procedures  are  documented  in 
[Ruch,  2001].  The  achieved  values  of  key  is 
via  transient  thermal  diffusivity  measurement 
[Appendix  A]. 

The  Table  shows  that  there  was  slight  wire 
diameter  growth,  caused  by  a  combination  of 
wire  stretching  during  the  weaving  and  metal 
addition  due  to  the  brazing  process.  More 
significant  is  the  decrease  in  weave  mesh 
numbers,  indicative  of  a  “looser”  mesh  than 
targeted.  This  is  due  to  the  inability  of  the 
weaver  to  maintain  a  “tight”  weave,  caused  by 
the  stiffness  of  the  metal  wire  filaments  and 
their  breakage  if  too  much  tension  is  applied  to  the  wire  filaments  during  the  weaving  process. 
The  loosening  of  the  mesh  results  in  an  increase  in  porosity,  and  a  decrease  in  the  number  of  wire 
intersections  that  can  be  brazed  together.  It  was  found  that  about  80%  of  wire  filament 
intersections  are  satisfactorily/partly  bonded.  The  remainder  are  either  un-bonded,  or  the  bonds 
contain  faults  [Ruch,  2001].  The  reduction  in  metal  fraction  and  wire  intersection  bonds  gives 
rise  to  a  reduction  in  heat  transfer  surface  area  (P),  and  effective  thermal  conductivity  (key),  and 
this  results  in  a  reduction  on  performance  (Eq.  10). 

Three-dimension  weave  setup  is  very  labor  intensive,  and  it  must  be  carefully  controlled  since  the 
weaving  of  aluminum  filaments  is  problematic.  Relatively  small  errors  in  mesh  number  result  in 
significant  changes  in  porosity,  specific  surface  area  and  effective  thermal  conductivity.  As  a 
consequence,  this  methodology  is  probably  best  applied  to  situations  where  large  volume 
manufacture  is  anticipated. 

Porosity,  Specific  Surface  Area  and  Effective  Thermal  Conductivity 

Figure  35  shows  a  unit  cell  of  the  weave  of  Fig.  1.  The  wire  pitch  in  each  of  the  coordinate 
directions  are  Mx~l ,  My~x  and  Mz'x  (M  is  the  “mesh  number”).  Consideration  of  this  unit-cell 
shows  that  the  porosity  (s)  and  specific  surface  area  (P)  are  given  by  the  following  expressions: 


Fig.  35  Three-D  stacked  weave  unit  cell. 
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(49) 


1  ~e-~7  ipx2MyMz  +  dy2MxMz  +  dz2MxMy) 

P  =  7t(dxMyMz+  dyMxMz+  dzMxMy)  (50) 

The  quantity  (1-s)  is  the  volume  “metal  fraction”. 

For  the  current  application,  we  require  that  wire  filaments  touch  at  all  possible  filament 
intersections  so  that  filaments  can  be  bonded  to  facilitate  conduction  within  the  mesh.  Then 


Mx  =  • 


1 


dy  +  dz 


My  = 


1 


dx  +  dz 


Mz-- 


dx  +  dy 


(51) 


Furthermore,  let  dx  -  dz,  and  define  the  filament  diameter  ratio,  r  =  dy/dx.  Equations  (49)  and 
(50)  become 


i  * 

l-e  =  — 
4 


1  +  r  +  r2 


(52) 


Pdy  =  S{\-e) 

where 


(53) 


S  = 


4  r 


l  +  2r 
1  +  r  +  r2 


(54) 


is  the  shape  factor.  Equations  (52  -  54)  show 
that  once  we  restrict  our  attention  to  bondable 
3-D  weaves;  the  yft/v-product  and  the  porosity 
are  functions  of  the  filament  diameter  ratio 
only. 

Figure  36  plots  fidy  of  the  3D  stacked  weave 
versus  metal  fraction  for  diameter  ratio,  0  <r 
<  <x>.  The  figure  shows  that  bondable  3-D 
stacked-weaves  with  r  >  1  can  be  configured 
to  have  metal  fraction  ranging  from  0.589  to 
jt/4  while  2.36  <  (5dy  <  2n.  Our  design 
objective  is  to  fabricate  a  3D  stacked  weave 
heat  exchange  matrix  having  r  =  2.  The 


Fig.  36  Specific  surface  area  vs.  metal 
fraction  of  porous  materials. 
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target  metal  fraction  is  then  (1-e)  =  0.61,  and 
P  =  4581m1.  The  target  point,  {r  -  2,  (1-e)  = 
0.61,  pdy  =  3.49}  is  also  shown  in  the  in  the 
figure. 

Furthermore,  the  specific  surface  area  formula, 
Eq.  (53)  takes  the  familiar  form  of  other 
porous  media.  For  example,  the  shape  factor 
for  an  unconsolidated  bed  of  spheres  is  S  =  6 
[Kaviany,  1995]  and  Xu  and  Wirtz  [2003] 
have  shown  that  isotropic  plain-weave  screen- 
laminates  (dy=dz,  My-Mz)  have  S  =  4.  These 
functions  are  shown  in  the  figure.  The  figure 
shows  that  our  target  design  point  will  have  a 
porosity  and  specific  surface  area  roughly 
equivalent  to  that  of  an  unconsolidated  bed  of 
spheres.  Larger  diameter  ratios,  r  will  result  in 


larger  achievable  specific  surface  areas. 


Xu  and  Wirtz  [2003]  show  that  isotropic  plain-weave  laminates  can  be  fabricated  to  have  0<  (1-e) 
<  0.534.  Under  these  conditions,  0  <  pdy  <  2.13.  The  figure  shows  that  addition  of  a  third 
filament  to  create  a  3D  orthogonal  weave  allows  for  the  structuring  of  a  mesh  having  a  larger 
metal  fraction  and  considerably  larger  specific  surface  area. 


3D  Woven-Mesh  Effective  Thermal  Conductivity.  Following  Chang  [1990],  we  transform  the  x- 
and  z-direction  wire  filaments  of  the  unit  cell  into  rectangular  cross  section  segments  shown  in 


Fig.  37.  Each  rectangular  wire  filament  has 
thickness  jid  1 4g  and  width  gd  so  that  the 

cross  section  area  of  each  filament  is  to!2  /  4. 
We  further  require  that  the  volume  of  the  unit 
cell  be  preserved  across  the  geometric 
transformation.  Then  g  =  n  1 4 . 

The  in-plane  effective  thermal  conductivity  in 
the  y-direction  may  be  determined  by 
considering  the  thermal  circuit  for  conduction 
in  the  y-direction  across  the  transformed  unit 


Fig  38  Thermal  circuit  for  y-direction 
conduction. 
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cell,  shown  in  Fig.  38.  Referring  to  the  transformed  unit  cell:  the  parallel  path  ( Rsl,Rfl ) 
represents  conduction  along  the  axis  of  the  central  dy-filament  and  surrounding  fluid;  the  series 
paths  (R2  =  Rf2+Rs2 )  and  (R3  =  Rs3+RJ3)  represent  conduction  across  dx-  and  dz-filaments  and 
fluid  regions  above  and  below  each  filament,  respectively;  and,  R4  represents  conduction  across 
the  dx-dz  filament  intersections. 


R\  =  - 


4(dx  +  dy)  •  dty~ 


jiks-4km 


\dy  dy)  V  [dy  dy) 


R2  = 


dx  •  kf  +  dz  ■  ks 


K  ■  kf  ■  g  -  dx  [dz(\-  g)+  dy\ 


(55) 


(56) 


R3  = 


dz-kj  +dz-ks 

ks-kf-g  dz-[dx{i-g)+dy] 


(57) 


R4  = 


cbc  +  dz 

ks -(g2- dx-dz) 


(58) 


where  ks  and  kf  are  the  thermal  conductivity  of  the  solid  and  “fluid”  phases,  respectively.  If  we 
define  the  effective  thermal  conductivity  in  the  y-direction  as 


then 


key  (dx  +  dy)(dz  +  cty)AT 
dx  +  dz 


(59) 


Ke,= 


^-•r2(l-^)  +  4^[r2  +  (l-g)2]  +  8^-r-(l-g)  +  4g2  4k,-g-(r  +  l-g) 


4(1  +  rf 


(kfi+l)(r  +  l)2 


(60) 


k&  k 

where  %e  =  — -  and  K  =  —  •  hi  most  cases  the  ratio  kfS  is  small;  so,  if  k/s &0  the  dimensionless 
y  K  ks 

effective  thermal  conductivity  becomes  a  function  of  diameter  ratio,  r  only, 

Kel  =  — |  — 7— — rr—  I  (61) 


1  *1 

f r2  +x/4' 

1  ^ 

1 

{ (r+if  ) 
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Figure  39  plots  Key(kfs= 0)  vs  the  metal 
fraction.  Also  shown  in  the  figure  is  the 
expected  dimensionless  thermal  conductivity 
for  a  bed  of  fused  spheres,  metal  foam  and  that 
for  a  plane  weave  screen  laminate.  The  figure 
shows  that  for  r  >  1,  the  dimensionless 
effective  thermal  conductivity  of  the  3-D 
weave  can  range  from  0.351  up  to  tc/4  while 
the  metal  fraction  ranges  from  0.589  up  to  jt/4. 
The  target  design  point,  {r  =  2,  ( 1-s)  =  0.61} 
gives  Key  —  0.481.  This  is  compared  to  a 
fused  bed  of  spheres,  which  is  expected  to 
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Fig.  39  Effective  thermal  conductivity  of 
exnanded  materials. 


have  -Ke/spheres)  «  0.15.  Furthermore,  incorporation  of  a  third  wire  filament  results  in  a 
dramatic  increase  in  the  effective  thermal  conductivity  relative  to  that  obtained  with  screen 
laminate  structures. 

Heat  Transfer,  Pressure  Drop  Correlation 

We  postulate  that  the  mesh  heat  transfer  coefficient,  h  is  functionally  related  to  fluid  and  flow 
properties  as  follows: 


Dimensional  analysis  gives 


h  =  f{p,G,jLi,cftkfdy,p,t,r,s) 


St  =  >(Re,Pr,r,s,-^) 


(62) 


(63) 


Cr D  A. 

where  Re  =  — -  is  the  mesh  Reynolds  number  with  n  =  _£  the  mesh  hydraulic  diameter.  In  a 
M  "  P 

similar  way,  we  postulate  that  the  pressure  drop  across  the  mesh,  AP,  is  functionally  related  to 
fluid  and  flow  properties  as  follows: 


A P  =  fn{p,  G,  p,  p,  t,  r,  e) 
Then  dimensional  analysis  gives 

f  =  jh(Re,r,£) 


(64) 


(65) 
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where  f  =  —  js  the  friction 

J  G  2  fit 

factor.  It  is  noted  that  Eqs.  (62)  and 
(64)  are  the  same  as  Eqs.  (28)  and  (30) 
for  isotropic  plain-weave  laminates  (r 
=  1).  Therefore,  we  anticipate  that 
Eqs.  (63)  and  (65)  will  take  the  same 
functional  form  as  their  correlations  or 
friction  factor  and  Stanton  number. 

Experiments  are  performed  to 
measure  the  pressure  drop  and  porous 
wall  effective  conductance.  Then,  the 
mesh  Stanton  number  is  determined 
from  Eq.  (8).  The  functional  forms  of  Fig.40  Modified  friction  factor  and  Colburn  j-factor. 
Eqs.  (63)  and  (65)  are  then 

determined.  Pressure  drop  and  heat  transfer  experiments  use  the  same  channel  flow  apparatus 
and  data  reduction  procedures  as  described  in  Appendix  B. 

Figure  40  compares  measured  friction  factor, 

(66) 


/  = 


A P 

!£!, 
2  P 


and  j-factor. 


j  =  St-  Pr 


2/3 


(67) 


for  the  present  3D  stacked  weave  (s  =  0.49,  t  =  6.35mm,  9.53mm  and  12.7mm)  with  correlations 
for  spheres  (s  =  0.39)  and  isotropic  plain-weave  screen-laminates.  Ninety-five  percent  (95%) 
confidence  level  error  bars  are  shown  on  the  figure.  The  present  data  lays  intermediate  to  these 
two  other  configurations.  The  friction  factor  is  correlated  as 


/  = 


1.078 +  — 

fl-e} 

ReJ 

l  £  J 

x-0.4 


(68) 


The  first  term  of  Eq.  (68)  represents  the  inertial  contribution  to  the  pressure  drop;  and,  the  second 
term  represents  the  viscous  contribution.  Since  the  present  experiments  consider  only  a  single 
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Table  5  Three-dimensional  orthogonal  mesh  characteristics. 


Characteristic 

3D  Weave 

Screen-Laminate 

External 

Characteristics 

t  =  3.175mm,  2 :H=  2t,  4t,  8t 
Aluminum,  Alloy  1100 

Coolant  is  Air  @  300K  (Pr  =  0.69) 

Internal 

Geometric 

dy  =  0.41  mm 
r  =  2,e  =  0.385 
“stacked” 

Md  =  0.577 
£=0.466 
dy  =  dz  =  0.35  mm 
“in-line”  laminate 

ySK1] 

8589, Eq.  (5) 

6073,  Eq.  (25) 

kev  [W/mK] 

84,  Eq.  (13) 

44  W/mK,  Eq.  (26) 

Eq.  (24) 

Eq-  (32) 

Eq.  (25) 

Eq.  (34) 

mesh  porosity,  we  have  used  Part  et  al’s  [2002]  term  for  the  second  factor  in  Eq.  (68).  Equation 
(68)  reproduces  the  data  that  generated  it  with  a  standard  error  of  4.4%.  In  a  similar  manor,  the  j- 
factor  is  correlated  with  a  power-law 


7  =  1.89  Re'044 


\-£ 
K  £  . 


v-0.2 


(69) 


Equation  (69)  reproduces  the  data  that  generated  it  with  a  standard  error  of  10.7%. 

Exchange  Matrix  Performance  Evaluation 

Three-dimensional  stacked  weaves  offer  considerable  design  flexibility.  Adjustment  of  the  wire 
diameter  ratio,  r  allows  for  control  of  the  structure’s  porosity,  heat  transfer  surface  area  to  volume 
ratio,  and  effective  thermal  conductivity.  However,  the  friction  factors  and  Stanton  numbers  of 
the  heat  exchange  matrix  are  comparable  to  those  of  other  heat  transfer  surfaces.  The  question 
that  must  be  addressed  is:  under  what  conditions  does  the  3D-weave  technology  offer  superior 
performance. 

Recognizing  that  we  have  shown  that  screen-laminates  can  generally  be  configured  to  out¬ 
perform  fused  particle  systems,  in  the  following  we  describe  a  fixed  outer  geometry  comparison 
of  the  performance  of  a  3D-weave  exchange  matrix  with  a  screen-laminate  exchange  matrix  of 
the  same  volume  and  face  area.  The  two  systems  are  deployed  as  in  a  single-fluid  parallel-plate 
heat  exchanger  such  as  a  cold-plate  or  flow-through  module. 

Table  5  summarizes  the  characteristics  of  the  two  systems.  Both  exchange  matrices  are  made  of 
aluminum  (alloy  1 100).  They  have  the  same  thickness,  face  area  and  plate-to-plate  spacing  (three 
spacing  are  considered).  The  coolant  is  air.  Consideration  of  higher  Prandtl  number  fluids 
produces  similar  results.  The  3D  weave  is  made  up  of  0.4 1  mm/0.20mm  ( 1 6mil/8mil)  wire  (r  =  2) 
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resulting  in  an  exchange  matrix  with  a 
porosity  of  0.385,  specific  surface 
area  of  8589  m'1  and  effective  thermal 
conductivity  of  84  W/mK. 

The  tightest  weave  achievable  with  an 
isotropic  plain-weave  has  the  product 
Md  =  0.577.  This  results  in  a  porosity 
of  0.466.  A  mesh  number  equivalent 
to  the  3d  weave  has  M  —  16.4  cm'1. 

Then,  the  wire  diameter  must  be  Fig.  41  Effective  conductance  and  superficial  coolant 
0.35mm  (0.014  in).  With  no  velocity  of  3D  stacked  weave.  Coolant  is  air  at  300K. 

interleaving  of  adjacent  screen  layers,  the  specific  surface  area  becomes  6073  m'1,  and  the 
effective  thermal  conductivity  is  44  W/mK,  roughly  half  the  value  achieved  with  the  3D  weave. 


Figure  41  plots  the  effective  surface  conductance  of  the  3D  stacked  weave  as  a  function  pressure 
drop.  Results  for  three  plate-to-plate  spacings  are  shown.  Also  shown  in  the  figure  is  the 
superficial  velocity,  V„  as  a  function  of  applied  pressure  drop.  At  AP  =  1250  Pa  (5  in  H20),  the 
conductance  exceeds  6500  W/m2K,  a  value  normally  associated  with  liquid-flow  turbulent 
convection  or  phase-change  heat  transfer.  At  this  pressure  drop,  the  superficial  coolant  velocity  is 
approximately  2.7  m/sec. 


We  define  the  pressure  drop  ratio  and  “heat 
duty”  ratio  as 


^  - 


AP(3D  •  weave) 
AP(screen  ■  la  min  ate) 


(26a) 


U  (3D  ■  weave) 

K  = - 

U ( screen  •  la  min  ate) 

(26b) 

Figure  42  plots  the  pressure  drop  ratio  and 
heat  duty  ratio  as  a  function  of  superficial 
coolant  velocity,  comparing  the  3D-weave  to 
the  screen  laminate.  In  eveiy  case,  the  screen 
laminate  pressure  drop  is  significantly  lower 
than  that  of  the  3D  weave.  However,  the 


10 


Rap 


. 1 2H  =  6.^5  mm 

-  2H=  12.7  mm 

- 2H  =  25.4  mm 


Fig.42  PEC  3D  stacked  weave  vs  isotropic 
plain-weave  screen-laminate. 
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corresponding  heat  duty  ratios  are  greater  than  one,  indicating  that  the  3D  weave  out  performs  the 
screen  laminate  system  by  about  50%.  Figure  42  shows  that  the  heat  duty  ratio  increases  with 
increasing  superficial  velocity.  This  is  a  consequence  of  the  3D-weave  having  a  higher  specific 
surface  area  and  mesh  heat  transfer  coefficient.  The  figure  also  shows  that  the  heat  duty  ratio  also 
increases  with  increasing  array  height.  This  is  a  consequence  of  the  significantly  higher  effective 
thermal  conductivity  achieved  with  the  3D-weave. 
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Conclusions 


A  mathematical  model  of  the  thermal  performance  of  a  woven  mesh,  deployed  as  a  heat  exchange 
surface,  is  developed.  Woven  meshes  are  modeled  as  porous  media.  Since  the  fluid  flow  path 
length  through  the  porous  wall  is  short  (by  design)  and  flow  rates  are  relatively  high,  a  two- 
energy  equation  model  is  developed.  The  present  model  shows  that  the  effective  conductance  of  a 
thin,  porous  heat  exchange  matrix  is  proportional  to  the  square  root  of  the  product  of  the  specific 
surface  area,  effective  thermal  conductivity  and  mesh  heat  transfer  coefficient.  Furthermore,  we 
propose  that  mesh  heat  exchange  matrices  be  laid  out  in  a  serpentine  pattern,  similar  to  the 
construction  of  a  conventional  automotive  air  filter.  This  configuration  maximizes  heat  transfer 
surface  area  while  it  minimizes  coolant  pressure  drop  (due  to  the  reduced  superficial  mass 
velocity).  Analysis  of  a  screen-fin  heat  sink  indicates  that  the  optimal  convergence/divergence 
angle  of  the  serpentine  layout  is  «  20°  -  30°. 

Screen  Laminate  Technology  Screen  laminates  are  simple  to  fabricate  structures.  They  can  be 
constructed  to  have  a  wide  range  of  porosity,  heat  transfer  surface  area  and  effective  thermal 
conductivity.  Isotropic  orthogonal-weave  laminations  can  be  configured  to  have  metal  fraction 
(l-«)  and  specific  surface  area  ranging  up  to  0.534c/-"1  and  2.13  ( cf  d)  respectively.  The  in- 
plane  effective  thermal  conductivity  can  be  as  much  as  6-  times  greater  than  the  cross-plane 
effective  thermal  conductivity,  and  the  in-plane  effective  thermal  conductivity  can  range  to  19.2% 
of  the  filament  material  value  (a  thirty  to  fifty  percent  improvement  over  other  porous  materials). 
Diamond  weaves  (and  anisotropic  weaves)  can  be  configured  to  have  proportionately  higher 
material  properties,  with  the  in-plane  effective  thermal  conductivity  ranging  to  78.5%  of  the 
filament  material  values. 

Dimensional  analysis  is  used  to  establish  the  definitions  of  the  dimensionless  groups  that 
characterize  the  pressure  drop  and  mesh  heat  transfer  coefficient.  Experiments  with  air  and  water 
have  been  used  to  establish  correlations  for  the  friction  factor  and  mesh  Stanton  number.  Screen 
laminates  are  found  to  exhibit  orthogonal  flow  pressure  drop  and  mesh  heat  transfer  coefficients 
comparable  to  other  porous  media.  Experiments  with  non-orthogonal  flow  (superficial  flow 
incidence  angle  =  (p)  show  that  the  friction  factor  increases  as  <p  -»  0;  and,  the  mesh  heat  transfer 
coefficient  increases  as  (p  — >  20°. 

A  performance  comparison  shows  that  screen-laminate  systems  can  generally  be  configured  to 
offer  thermal  and  pressure  drop  performance  superior  to  unconsolidated  packed  bed  matrices. 
This  is  particularly  true  for  exchange  matrix  configurations  that  exploit  the  higher  effective 
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thermal  conductivity  achievable  with  screen-laminate  systems.  A  design  exercise  of  a  screen-fin 
heat  sink  for  a  P4  microprocessor  cooling  application  shows  that  a  simple  device  can  be 
configured,  which  rivals  die  performance  of  a  dense-fin  (expensive)  unit  that  is  the  unofficial 
industry  benchmark. 

Stacked  Weave  Technology  Three-dimensional  stacked  weaves  are  relatively  dense  structures, 
with  the  y%/y-product  and  porosity  functions  of  the  wire  diameter  ratio  only.  Metal  fractions  can 
range  from  0.589  up  to  0.785;  and,  the  pdy- product  can  range  up  to  2n.  Metallic  weaves  can  be 
structured  to  have  effective  thermal  conductivity  that  is  two,  or  more  times  greater  than  what  can 
be  achieved  with  other  porous  media,  with  y-direction  effective  thermal  conductivity  approaching 
78%  of  the  filament  material  value.  Mesh  heat  transfer  coefficients  and  friction  factors  are 
comparable  to  those  achieved  with  other  expanded  materials.  However,  high  P-values,  coupled 
with  high  effective  thermal  conductivity  result  in  exchange  matrices  that  out-perform  other 
exchange  matrix  configurations  in  applications  where  coolant  pressure  drop  is  not  a  critical 
constraint. 

Three-dimension  weave  setup  is  very  labor  intensive,  and  and  it  must  be  carefully  controlled 
since  the  weaving  of  aluminum  filaments  is  problematic.  Relatively  small  errors  in  mesh  number 
result  in  significant  changes  in  porosity,  specific  surface  area  and  effective  thermal  conductivity. 
As  a  consequence,  this  methodology  is  probably  best  applied  to  situations  where  large  volume 
manufacture  is  anticipated. 
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Appendix  A:  Measurement  of  ke 


Consider  a  uniform  cross  section,  long  and  slender  test  article 
with  well-insulated  periphery,  as  shown  in  Fig.  Al.  The  test 
article  is  initially  at  a  uniform  temperature,  T*  when  a  time- 
varying  heat  flux  is  applied  to  one  end.  We  measure  the 
temperature  response  at  three  locations,  as  shown  in  the  figure. 
Assume  constant  thermophysical  properties  and  one¬ 
dimensional  transient  conduction  in  the  domain  between  points 
1  and  3.  The  temperature  response  at  point  2  will  be  given  as 

T2(t)  =  fct[T\(T),T3(T\AyhAy2,a)  (18) 


The  temperature  response  at  point  1  describes  the  heat  input  to  Fig.  Al  Test  Configuration, 
the  test  domain;  that  at  point  3  describes  the  heat  outflow  from  the  test  domain.  Then,  the 
temperature  response  at  point  2  can  be  used  to  determine  the  thermal  diffusivity,  a,  [7]. 


The  test  rig  consists  of  a  radiant  heat  source 
and  shutter  mechanism,  and  well-insulated 
screen  laminate  test  articles  that  are 
approximately  200mm  long,  as  shown  in  Fig 
A2.  Temperatures  at  points  1,  2  and  3  are 
sensed  with  0.25mm  diameter  Type  T 
thermocouples  that  are  soldered  to  the  test 
article.  Measurement  points  are  located  with 
an  accuracy  of  ±0.1  mm.  Temperatures  are 
recorded  at  a  rate  of  10  -  20  samples  per 
second,  with  an  accuracy  of  approximately 


Heat  Receiving 
Plate 

Foam 


strumented 
Test  Article 

f ermocouple 


Fig  A2  Insulated  test  articles. 


±0.2°C. 

Measurement  Calibration  Three  bar-shaped 
oxygen-free  copper  test  articles  (Alloy 
C10100,  a  =  ll.lxl0'5m2/s)  were  tested. 
Each  test  was  repeated  six  times.  The  mean 
measured  thermal  diffusivity  for  each  of  the 
test  articles  was  11.76xl0'5m2/s,  11.62x10' 
snt2/s  and  1 1.72x1  (T^rrf/s,  giving  an  overall 


Fig  A3  Copper  calibration  test  articles. 
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average  value  of  11.7xl0'5m2/s.  The  95%  confidence  interval  for  all  18  measurements  (2a  about 
mean  values)  is  0.7x1 0‘5m2/s  and  the  implied  offset  error,  based  on  the  handbook  value  of  thermal 
diffusivity  of  alloy  C 1 0 1 00,  is  +5.2%. 


Screen  Laminates  Test  Articles  Isotropic, 
orthogonal  weave  and  diamond  weave  plain- 
weave  copper  screens  are  stacked  and  bonded 
together  (using  95/5  Sn/Pb  solder)  to  form  the 
screen  laminates.  Various  screen  laminate 
samples  were  fabricated  with  this  method:  the 
mesh  numbers  range  from  6.3cm1  (Minch1)  to 
15.75cm1  (40inch1),  the  bare  copper  wire 
filament  diameter  ranges  from  0.28mm(l  1/100 
inch)  to  0.46mm(18/1000inch),  and  the 
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Fig  A4  Typical  test  result. 


number  of  layers  ranges  from  3  to  10.  Measured  compression  factors  range  from  0.70  to  1.00. 


Figure  A4  shows  a  typical  experimental  result  where  Tc(t;a)  is  the  calculated  temperature 
response  of  point  2  with  ae  =  5.35x1  O'5  m2/sec. 


In  the  calculation  of  ke,  the  effective  density/specific  heat  product  of  the  screen  laminate,  pc  is 
calculated  based  on  the  measured  weight  ratios  of  copper  and  solder  material  in  the  composites. 
The  copper  wire  properties  are:  p  =  8940  kg/m3,  c  =  393.5  J/(kg-K),  ks  =  400  W/(mK).  The  Sn-Pb 
solder  material’s  properties  are:  p  =  73 17kg/m3,  c  =  230.0J(/kg-K),  lc,  =  56  W/(mK).  As  different 
specimens  have  different  weight  ratios  of  copper  and  solder  material,  so  the  calculated  effective 
densities  for  the  laminates  range  from  8694kg/m3  to  8796kg/m3,  the  specific  heats  range  from 


369.1J/kg  K  to  379.0J/kg  K  and,  the 
thermal  conductivity  range  from 
349.7  W/mK  to  371.1  W/mK.  The 
thermal  conductivity  of  air 
(0.0278W/(m-K))  at  temperature  of 
323K,  is  used. 

Test  results.  The  average  test  result 
of  dimensionless  reduced  in-plane 
effective  conductivity  for 
orthogonal  weave  specimens  is 


Cf-(l-E) 


Fig.  A5  Benchmark  measurements  of  key  in  air. 
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shown  with  2 a  error  bars  in  Figure  A5.  Each  data  point  shown  represents  six  or  more  test  runs 
of  the  same  sample.  The  model  prediction  for  the  isotropic  plain-weave  screen  laminate  saturated 
with  air  is  also  shown  in  the  figure.  It  shows  that  the  prediction  of  the  present  model  is  within 
10%  of  all  of  the  measurements.  With  consideration  of  the  non-uniform  solder  layer’s  thickness 
on  the  filaments,  and  non-uniform  properties  of  the  high  porosity  screens,  the  prediction  of  the 
present  model  is  quite  accurate. 
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Appendix  B:  Heat  Transfer/Pressure  Drop  Measurements 


Pressure  drop  and  heat  transfer  experiments 
are  done  in  two  different  channel-flow 
apparatus.  A  schematic  of  the  test  section  for 
either  apparatus  is  shown  in  Figure  Bl.  The 
figure  shows  an  edge  view  of  a  screen 
laminate,  of  thickness  t,  located  in  a  parallel- 
plate  channel,  which  is  approximately  18  mm 
high  x  100  mm  wide.  A  fluid,  at  mass  velocity 
(GO  and  temperature  (Tf>j)  passes  through  the 
test  article.  The  test  article  is  heated  from  the 
side.  Thermocouples  measure  the  inlet 
temperature  and  base  temperature.  Coolant 
pressure  drop  is  measured  with  an  electronic 
manometer  having  ±4%  accuracy. 


Fig.  Bl  Edge  view  of  mesh  exchange  matrix 
test  section. 


Air flow  experiments  The  channel  is  of  open- 
loop,  induced-draft  design.  Laboratory  air 
passes  through  a  honeycomb  flow 
straightener;  the  test  article;  a  second  flow 
straightener;  a  plenum  chamber  and  suitably 
long  pipe  to  a  laminar  flow  element,  which 
measures  the  volumetric  flow  rate;  and,  then 
to  a  variable  speed  exhauster.  The  laminar 
flow  element  has  ±3%  accuracy. 

A  thin  film  heater  applies  power  to  the  bottom 
(one  side)  of  the  test  article.  Thermocouples 
are  used  to  monitor  the  base  temperature  of  the  test  article  and  inlet  air  temperature.  We  estimate 
that  heat  losses  out  of  the  bottom  of  the  test  rig  are  less  than  2%  of  the  input  power  (q), 
temperature  measurements  are  accurate  to  ±0.05  °C,  and  q  is  measured  to  ±4%. 


Liquid  flow  experiments .  Experiments  are  conducted  in  a  closed  loop  chilled  water  flow 
apparatus  shown  schematically  in  Fig.  B3.  The  test  rig  consists  of  a  pump/fluid  reservoir,  heat 
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copper  screen  with  d  =  0.48  mm  (0.019”),  M  =  7.9  cm'1  (20m'1)  with  s  =  0.677 .  A  total  of  1 8 
samples  were  tested  where  the  flow  rate  was  varied  over  the  following  range:  1.76  <,  Re  <,  300. 
All  experiments  were  performed  at  least  twice  (on  different  days),  to  check  the  repeatability  of 
the  data.  A  Monte  Carlo  error  propagation  simulation  indicates  the  following  95%-confidence 
level  tolerance  on  computed  results:  Re  less  than  +5%;  f  less  than  ±  1 5%;  St  less  than  ±10%. 
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